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We present upper limits on line emission in the Cosmic X-ray background (CXB) that would 
be produced by decay of sterile neutrino dark matter. We employ the spectra of the unresolved 
component of the CXB in the Chandra Deep Fields North and South obtained with the Chandra 
CCD detector in the E = 0.8 — 9 keV band. The expected decay flux comes from the dark matter 
on the lines of sight through the Milky Way galactic halo. Our constraints on the sterile neutrino 
decay rate are sensitive to the modeling of the Milky Way halo. The highest halo mass estimates 
provide a limit on the sterile neutrino mass of < 2.9 keV in the Dodelson-Widrow production 
model, while the lowest halo mass estimates provide the conservative limit of nis < 5.7 keV (2cr). 
We also discuss constraints from a short observation of the softer {E < 1 keV) X-ray background 
with a rocket-borne calorimeter by McCammon and collaborators. 



PACS numbers: 95.35.-|-d,14.60.Pq,14.60.St,98.65.-r 



I. INTRODUCTION 

The abundance of cosmological dark matter is now well 
quantified by cosmological observations to better than 
10% in density [HQj yet its identity remains unknown. 
Several indicators point towards a modification of the 
properties of dark matter clustering from the cold dark 
matter (CDM) paradigm to resolve potential discrepan- 
cies that persist with the ansatz of an absolutely cold 
dark matter candidate. Among potential problems with 
CDM are indications for density cores in local group 
dwarf spheroidal galaxies d, H, H, Q and the lack of a 
correspondence of the observed dwarf galaxies with the 
number of halos expected with the same maximal ve- 
locity dispersion or mass 0, Q • This has prompted the 
investigation of alternatives to CDM, such as warm dark 
matter (WDM) [1, late-decaying super WIMP dark 
matte r [111 . [T^ . fuzzy cold dark matter [l^, or meta- 
CDM [ij^ 

A leading particle candidate for WDM is a fermion 
that has no standard model interactions, yet couples with 
the standard model neutrino sector via the neutrino mass 
generation mechanism, namely a sterile neutrino. Exten- 
sions to the standard model of particle physics typically 
include sterile neutrinos, including left-right symmetric 
(mirror) models [13], supersymmetric axinos as sterile 
neutrinos [l6| , superstring models [it'I , models with large 
extra dimensions llSL | l9ll . and phenomenological models 
such as the i^MSM [23|. Other motivations for a ster- 
ile neutrino to have parameters in the parameter space 
of interest for their oscillation production as dark matter 
are the generation of pulsar kicks 0, [23. l23j , Type II su- 
pernova shock heating enhancement [24!|, and enhanced 



molecular hydrogen formation at high redshift [25| . 

Sterile neutrinos may be produced in the early universe 
via non-resonant Dodelson-Widrow (DW) oscillation pro- 
duction [26| . or via a resonant oscillation production in 
cosmologies with a non-zero lepton number Both 
non-resonant DW and resonant production models fall 
in a parameter space that is continuous in the cosmolog- 
ical lepton number, with non-negligible lepton numbers 
allowing for resonant production at smaller mixing an- 
gles. The DW model is the simplest model for sterile 
neutrino dark matter production, because it assumes a 
standard thermal history in the early universe, zero lep- 
ton numbers, and no additional couplings of the sterile 
neutrino. For the non-resonant and resonant cases, pro- 
duction occurs near the quark-hadron transition, and the 
relation between the critical density in sterile neutrinos 
and its mixing parameters with active neutrinos is given 
by Refs. [1^, [2^, [s^l ■ Sterile neutrinos can be produced at 
smaller mixing angles than the DW mechanism through 
the resonant production model or via additional cou- 
plings to other particles, such as the inflaton 'sT]. They 
could be produced at larger mixing angles than the DW 
model if over-abundance due to the nonresonant oscilla- 
tion production is avoided through dilution by massive 
particle decay after production [32] or by a low reheating 
temperature scale in the early universe (33| . 

A particularly interesting particle mass range for the 
sterile neutrino can be framed by the mass required to 
produce a constant density core of 100-300 pc like that 
discussed for the Fornax dwarf spheroidal [1, which 
is 0.5-1.3 keV [3, [S^], or more massive if the phase- 
space packing limit is not achieved. In addition, the re- 
quirement for a viable sterile neutrino dark matter can- 
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didate constrains their parameter space by demanding 
that: 1) the total decay time-scale is larger than the 
age of the universe, 2) it is consistent with radiation en- 
ergy density constraints from primordial nucleosynthesis 
and the cosmic microwave background, 3) it is consis- 
tent with lithium photoproduction constraints, and 4) 
it is consistent with the diffuse X-ray background, and 
X-ray observations of of clusters or other nearby struc- 
tures [H, m, [131 ■ Sterile neutrinos are constrained by 
X-ray observations because of the considerable radiative 
decay width to a photon and lighter-mass neutrino due 
to the same coupling required for their production, pro- 
ducing a spectral line in the X-ray [111 . 

Abazajian, Fuller & Tucker [s^ found that X-ray ob- 
servations could present the most stringent constraints in 
the upper particle-mass range regions of parameter space. 
An estimate in that work of the sensitivities of X-ray ob- 
servations of the large dark matter overdensities present 
in clusters of galaxies and field galaxies was made and led 
to an estimated limit from XMM-Newton observations 
of the Virgo cluster on the non-resonant DW production 
model on the particle mass of the sterile neutrino. Recent 
work has shown that the local radiative decay flux from 
dark matter overdensities in Local Group structures can 
be comparable to that estimated from clusters of galax- 
ies and field galaxies, with reduced continuum emission 
and bigger angular size, and therefore increased sensitiv- 
ity [39|. This was also found in subsequent observational 
work [4dl [4l[ . [i^ . For example, using a nondetection 
of the decay line from the Andromeda galaxy by XMM- 
Newton, Watson et al. [iO] derived vris < 3.5 keV in the 
DW model (95% CL). 

Upper mass constraints from the lack of X-ray line flux 
are complementary to lower mass constraints from ob- 
served cosmological structure. Since sterile neutrinos be- 
have as increasingly warm dark matter for lighter particle 
masses, observations of the lack of deviations from abso- 
lute cold dark matter clustering constrain lighter parti- 
cle masses. Two recent analyses of a single measurement 
of the flux power spectrum of the Lya forest from the 
Sloan Digital Sky Survey (SDSS) [H find roughlysimi- 
lar constraints rn^ > 14 keV [31 and > 9 keV ^ in 
the DW model. These improved considerably relative to 
previous constraints at rris > 2 keV, which used the in- 
ferred linear matter power spectrum from the SDSS and 
higher-resolution Lya flux power spectra [i^ |43|. The 
improvement by a factor of five stems from the high- 
redshift (z ~ 4) Lya flux power spectra of the SDSS, 
where there is less enhancement of the amplitude of mat- 
ter power at small scales due to the nonlinear growth 
of structure. The newer Lya limits are quite stringent. 
When combined with X-ray observations of Andromeda 



^ We have applied a rough 10% reduction in the limit of Ref. l45j 
in order to include relevant nonthermal effects, though this cor- 
rection is significantly mass dependent [Ssl . 



and other mass halos, they would exclude the DW pro- 
duction mechanism for sterile neutrinos. However, it 
should be noted that the amplitude and slope of the dark 
matter power spectrum inferred from the flux power spec- 
trum measurement of Ref. [43| used in both sterile neu- 
trino analyses is inconsistent with that inferred from the 
Wilkinson Microwave Anisotropy Probe (WMAP) third 
year data [l| and also indicates a number of relativistic 
degrees of freedom of N^, = 5.4lg;g which is in ten- 
sion both with that expected N^, = 3.046 [4§| for the 
case of active neutrinos alone and that constrained by 
primordial nucleosynthesis, = S.OSlto^gs This 
may indicate hidden systematic effects within the mea- 
surement presented in Ref. [i^ of the Lya flux power 
spectrum which could alleviate or remove constraints on 
the inferred matter power spectrum from the Lya for- 
est [53. 

In this paper, we use Chandra spectra of the unresolved 
component of the Cosmic X-ray Background (CXB) in 
the direction of Chandra Deep Fields North and South 
(hereafter CDFN and CDFS) [H [H to search for the 
contribution from the radiative decay of a sterile neutrino 
which could comprise the dark matter halo of the Milky 
Way (MW). In addition, we consider the sensitivity of 
an existing short measurement at lower energies (0.4-1 
keV) with an X-ray calorimeter [53| to the sterile neutrino 
signal. 



II. THE DARK MATTER MODEL 

For a sterile neutrino of mass rus and a mixing an- 
gle 9, the decay rate for a Dirac-type active-sterile mass 
coupling is given by [H, [ssj 



1.36x10"^^ s"^ 
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Note that here we identify the particle mass nig with that 
of the mass eigenstate most closely associated with the 
sterile neutrino. The dominantly sterile neutrino mass 
eigenstate decays into a photon of E — ms/2 and a pre- 
dominantly active neutrino mass eigenstate. The decay 
rate as a function of particle mass, for a ratio of dark 
matter density density to critical density Qg and a fixed 
quark-hadron transition temperature is given by 



r°W(ms) = 9.95x10 



-30 g-1 



1 keV 



3.37 



0.26 



(2) 



Here, we use the inferred relation between mass and mix- 
ing angle which is necessary to produce the appropriate 
cosmological density of dark matter in sterile neutrinos 
through the DW mechanism 28|. We assume a cross- 
over transition for the quark-hadron at a temperature of 
170 MeV for the relatoin, Eq. In general, this can 
vary significantly due to uncertainties in the nature of the 
quark-hadron transition, as pointed out in Abazajian & 
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Fuller '29'| and shown in greater detail recently by Asaka, 
Laine & Shaposhnikov ^561] . 

The expected flux from the decay of sterile neutrino 
dark matter depends on the mass of the particle, its decay 
rate, as well as the distance and distribution of dark mat- 
ter across the field of view of a detector, such as Chan- 
dra. Estimates of the dark matter mass within the field 
of view of the detector come from dynamic measures of 
the galactic mass profiles, or for clusters of galaxies, from 
the X-ray hydrostatic method or gravitational lensing. 

The flux in the neutrino decay line reaching the detec- 
tor is calculated as follows. Denoting x the linear coor- 
dinate along the line of sight and dfl an element of the 
detector field of view (FOV) , each corresponding volume 
element dfl dx contains pd-ai{x) /nis neutrinos, each de- 
caying with a frequency given by Eqs. (HH^]). In Eu- 
clidean geometry (i.e., as long as the object redshift z 
satisfies (1 -I- z)^ « 1), the photon flux from this volume 
element reaching a unit effective area of the detector is 



df ~ X dQ dx 



(3) 



The flux from the cone subtended by the detector FOV 
(in photons s""'^ cm~^) is therefore 



/ 



Aums 



dil / pdm(a^) dx 



dm? 



(4) 



where the first integral is over the FOV and the second 
integral is along the line of sight, and 5dm is the dark 
matter mass column density averaged over the FOV Q. 
The fiux coming from each direction in the FOV depends 
only on Sdm m that direction; details of the line of sight 
mass distribution, and even the distance to the object, 
are not important for objects at low z. For Chandra^s 
small FOV (a r = 5' circle used in this work), we can 
ignore the changes of the Milky Way mass column density 
across the FOV. However, in HIV] we will consider an 
instrument with a much wider FOV, for which we directly 
integrate the column density within the FOV. 

To calculate Sdm for the Milky Way halo in the di- 
rections of the two Chandra CXB observations, we use 
mass models from Galactic dynamics measures presented 
in Klypin et al. [s^l and Battaglia et al. [SSj, which are 
both consistent with each other and with Re f. [59| . We 
use the Navarro-Frenk- White (NFW) profile 60] 



p(r) 



1 



(5) 



where ps is the characteristic density and is the scale 
radius. The NFW profile is expected for CDM and should 
be consistent with WDM models where cores are ex- 
pected to be much smaller than the MW NFW scale 
radius [6l|. 

The dynamical constraints on the dark matter halo of 
the MW are complicated by the fact that we sit within 
this galaxy, and depend strongly on the measured or- 
bital velocities of satellite dwarf galaxies. The 68% con- 
fidence region of the NFW models from Battaglia et 



al. also covers the allowed model range from Klypin et 
al., while still being consistent with observations. We 
thus adopt that interval for our estimates. These mod- 
els have Mvir = 0.6 x 10^^ Mq (low mass MW) and 
Mvir = 2.0 X 10^2 (ijigi^ jj^ass MW), i?vir = 255 kpc, 
Rvir/fs = 18, and Rq — 8 kpc. In both models, most of 
the baryonic mass is concentrated in the central region 
of the MW (disk and bulge), which our lines of sight 
miss, thus we can simply (and somewhat conservatively) 
multiply the total density with the universal dark mat- 
ter fraction /dm — f^DM/(f^DM + ^b) ~ 0.867, where the 
fraction of critical density of the dark matter we take is 
f^DM = 0.26, and baryon density Of, = 0.04. The un- 
certainties in the MW halo model are much larger than 
that of the universal dark matter fraction, which we keep 
fixed. 

For a point at a distance x from the Sun along the 
line of sight with Galactic coordinates (£, h) , the distance 
from the center of the halo is 



(x^ - 2xRq cosb cos £ + R%y 



(6) 



which is used to evaluate the column density integral in 
Eq. g]). In the direction of the CDFN at £ = 125.89°, 
b = 54.83°, the dark matter surface density is 



0.0362 g cm-2 (high mass MW) 
0.0109 gcm-2 (low mass MW). 



For the CDFS at ^ = 223.57°, b = -54.44°, the corre- 
sponding surface densities are 4% lower. 

We integrated the MW halo out to i?vir. However, 
most of the column density accumulates at small radii 
(within O.liJvir), so the outer radial cutoff does not mat- 
ter. It also means that the column density depends on 
small-scale deviations from the symmetric model. There- 
fore, our range of dark matter column densities is only 
qualitative and does not reflect the full uncertainty, which 
is difficult to quantify, and therefore we leave our in- 
ferred limits with as they are but with the caveat of po- 
tentially non-negligible systematic uncertainty. A sub- 
stantive measure of the systematic uncertainty may be 
estimated via comparisons with other independent mea- 
surements. 

As a comparison, the mass surface density within 5° 
of the Galactic center is an order of magnitude higher, 
as is the surface density on the line of sight near the 
center of M31. Massive galaxy clusters have column den- 
sities of 0.1 — 0.3 gcm~^ in their core regions. Clearly, 
our lines of sight are not optimal for the neutrino line 
search, and more sensitive limits can be obtained, for ex- 
ample, from observations of the Galactic Center (which 
we will attempt in a future work). At the same time, 
our CDF datasets have the advantage of an almost com- 
plete removal of the sources causing the CXB and a very 
accurate background modeling. 
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FIG. 1: The three unresolved CXB spectra [CDFN-VF 
(black), CDFN-F (red), and CDFS (blue)] in the 0.4-10 keV 
band. They are fit well by a simple model consisting of a 
power law absorbed by the Galactic hydrogen column repre- 
senting the extragalactic component, plus unabsorbed local 
warm thermal emission. All three spectra are fit well, with 
a total reduced ^ of 0.87, and no residual line-like features 
are seen. Only energies E > 0.8 keV are used for the line flux 
limits. See text for details. 



III. CHANDRA CXB SPECTRUM 

The Chandra X-ray Observatory [6^ has a 1 arcsecond 
on-axis angular resolution and is uniquely suited for re- 
solving the CXB covered by its ACIS detector, which op- 
erates in the 0.4-10 keV energy band. It has performed 
a series of very deep observations of two fields, CDFN 
and CDFS, for a total exposure of 2 Ms and 1 Ms, re- 
spectively, aimed at resolving as much of the extragalac- 
tic CXB into point sources as possible [63l H^l • Hickox 
& Markevitch 52] (hereafter H06) have used these ob- 
servations, along with recent accurate calibration of the 
ACIS internal background, to derive a spectrum of the 
CXB that is still unresolved after spatially excising all 
X-ray sources detectable in these deep exposures. They 
found that about 20% of the CXB in the 1-2 and 2-8 
keV bands remains unresolved, while at £^ < 1 keV, the 
Galaxy contributes a dominant genuinely diffuse compo- 
nent. The 1-8 keV unresolved spectrum is well modeled 
by a power law with a photon index T w 1.5, and the 
Galactic diffuse component by thermal emission from a 
T « 0.2 keV plasma. 

These unresolved spectra do not exhibit any emission 
lines that can be attributed to decaying sterile neutrinos 
in the Galactic halo, and we will use this fact to place 
upper limits on the neutrino line in the E 0.8 — 8 keV 
range. We note that in a later work, Hickox & Marke- 
vitch ^] have additionally excised sources not detected 
in X-rays but seen by the HubbleSpa.ce Telescope in the 
optical. These sources collectively account for most of the 
unresolved 1-8 keV flux, with a residual CXB brightness 
consistent with zero. Although this seems advantageous 
for our upper limits, we will not use these new spectra. 



because excising those optical sources reduced the solid 
angle from which the CXB spectra are collected by a 
factor of 3, compared to that in HOG. This increased 
statistical uncertainty and removed any advantage for 
our analysis. The CXB in either case is a small frac- 
tion of the detector background, and our constraints on 
the neutrino line are limited by photon statistics of this 
background. We thus elected to tolerate some unresolved 
CXB flux (modeled with a power law as described below) 
but have a significantly smaller statistical scatter in each 
spectral bin. At the same time, we will take advantage of 
a longer calibration exposure of the ACIS internal back- 
ground that was used in the latter work. 

Full account of the CDF spectra derivation is given 
in HOG, and here we give only the relevant details. The 
CFDN dataset is divided into two, one observed with 
ACIS in FAINT mode and another in VFAINT (the latter 
more advantageous for background modeling). We treat 
these subsets as independent observations (CDFN-F and 
CDFN-VF). After excluding periods of elevated detec- 
tor background, their exposures are 472 ks and 537 ks, 
respectively. The CDFS was observed in FAINT mode 
and has a clean exposure of 5G8 ks. Solid angles, after 
the source exclusion, subtended by the CDFN and CDFS 
fields are 0.0135 deg^ and 0.0159 deg^, respectively. The 
internal detector background is modeled and subtracted 
using a set of calibration observations in which ACIS was 
shielded from the sky. We will use a more recent 325 ks 
calibration dataset from Hickox & Markevitch [5^ , com- 
pared to 23G ks in HOG. This exposure is still shorter than 
any of the CDF exposures, and the statistical uncertainty 
of this calibration spectrum will be the main limiting fac- 
tor for our analysis. Otherwise, the CXB spectra that we 
analyze here are identical to those in HOG. 

As noted in HOG, the three unresolved CXB spectra 
(CDFN-VF, CDFN-F, and CDFS) in the 0.4-10 keV 
band are fit well by a simple model consisting of a power 
law absorbed by the Galactic hydrogen column {Nh — 
1.5xl02"cm2 for CDFN and 0.9x10^^^ for CDFS) rep- 
resenting the extragalactic component, plus unabsorbed 
local warm thermal emission with solar heavy element 
abundances (the APEC model [11]). We assumed the 
power law to be the same between the CDFN-VF and 
CDFN-F spectra, but allowed it to be different for CDFS, 
since it points to a different region of the sky. Ther- 
mal models were allowed to be different between all three 
datasets, because they may also include a time-variable 
near-Earth charge exchange contribution (the dominant 
features of both the thermal and charge-exchange compo- 
nents is an Ovii line around 570 eV, so they are difficult 
to disentangle). 

The resulting fits are shown in Fig.[T] All three spectra 
are fit well, with a total reduced of 0.87; no residual 
line-like features are seen. We can place an upper limit 
on such a line as a function of energy, excluding the range 
-E < 0.8 keV, where such a limit would not be interesting 
because of strong emission lines in the Galactic thermal 
spectrum, not resolved well by the ACIS CCD. To do this. 
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FIG. 2; Shown are the hmits on the flux in a hne as a func- 
tion of energy from the CDFN-VF, CDFN-F, and CDFS. For 
each line energy in the range 0.8-9 keV, all model parameters 
were allowed to vary while deriving upper limits on the sterile 
neutrino line normalization. The upper line and (cyan) band 
are the 3a limit, with the band representing uncertainty in 
the detector background modeling. The lower band and line 
correspond to that from the 2a limits. The diagonal band is 
the range of expected line flux from the MW halo as a func- 
tion of energy for the sterile neutrino in a DW production 
model. 



we added a monochromatic line to the spectral model 
with the relative normalizations between the CDF re- 
gions determined from the predicted neutrino brightness 
in the directions of those fields ( pi)) and the region solid 
angles. For both the low-mass and high-mass Galaxy 
models, this ratio for CDFS/CDFN is 1.14. For each Une 
energy in the range 0.8-9 keV, we re-fit the model and 
derived upper limits on the neutrino line normalization, 
allowing all other model parameters to vary as above (so 
that the neutrino line was allowed to account for the flux 
from other model components). 

The spectra were binned in such a way that an emission 
line would be resolved at all energies and the statistics 
in each bin is Gaussian. The upper limits on the line 
normalization can then be derived from Ax^ w.r.t. the 
best- fit model. It is an adequate estimate (cf. Ref. [11]), 
because the sky signal is a small fraction of the underly- 
ing detector background count rate (20% and 3% in the 
1-2 and 2-8 keV bands, respectively), which dominates 
the statistics. We note that while the CDF spectra are 
statistically independent, the detector background spec- 
trum is almost the same for all three datasets, which 
introduces correlations among the three resulting CXB 
spectra. This has been taken into account by a simple 
Monte-Carlo simulation, from which we determined that 
for our particular combination of the CDF and back- 
ground exposures, 2a and 3cr upper limits for one in- 
teresting parameter correspond to Ax^ of 8.8 and 19.7, 



respectively (compared to 4 and 9 for uncorrelated spec- 
tra). Using these values, we obtained upper limits on the 
neutrino line flux shown in Fig. [2l The constraints gen- 
erally follow the energy dependence of the ACIS effective 
area and the presence of bright background lines (e.g., 
the fluorescent Au line at E = 2.1 keV). In addition, 
we varied the normalization of the detector background 
spectrum by ±3% (simultaneously for all three spectra), 
which represents the background modeling uncertainty 
(H06). The effect of this is shown as bands around the 
2 and 3cr limits in Fig. [2l We use the upper envelopes 
of these bands for the constraints below. The parameter 
space excluded by these constraints are shown, relative 
to other constraints, in Fig. [3l 

In the DW model, these upper limits correspond to 2a 
limits on the sterile neutrino particle mass 

f2.87keV (high mass MW) 
™^ ^ 1 5.66 keV (low mass MW). ^ ' 

As is clear from these limits, there is considerable uncer- 
tainty due to the modeling of the dynamics of the MW 
halo. 



IV. OBSERVATIONS OF THE SOFT X-RAY 
BACKGROUND 

Here we explore the constraints from a related obser- 
vation at £^ < 1 keV using a detector with a much higher 
spectral resolution. The observation was made with an 
X-ray calorimeter flown on a sounding rocket, as reported 
by McCammon et al. [53 |. The spectrum from their 
100 s exposure is shown in the upper panel of Fig. U) 
The observation was made towards Galactic coordinates 
£ — 90°, 5 — -f60°. The mass surface density averaged 
over the wide {il = 0.81 sr) field of view of this experi- 
ment is 

- [0.0460 g cm-2 (high mass MW) 

~ 1 0.0138 g cm-2 (low mass MW). ^ ' 

The quantity iliSdm, which determines the expected ster- 
ile neutrino decay signal in the spectrum, is 2.5 x 10^ 
times higher than in the CDF observations. The dwarf 
spheroidal galaxies Ursa Minor and Draco lie within the 
FOV of this observation, but contribute only 2% to 5dm- 
To get a rough idea of how sensitive this experiment is 
to the sterile neutrino decay, we estimated the neutrino 
line flux at E = 1 keV. Using the detector response and 
exposure time from McCammon et al., we find that an 
emission line from the nig = 2 keV sterile neutrino in the 
DW model would contain 1.8 or 0.5 counts in our maxi- 
mum and minimum halo mass models, respectively. This 
is obviously too low to be detectable in the present data. 
A rough upper limit on the line flux at E = 1 keV can be 
derived in the following way. The observed background 
in the McCammon spectrum at 1 keV is ~ 1.5 counts 
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FIG. 3: Full parameter space constraints for the sterile neu- 
trino production models, assuming sterile neutrinos consti- 
tute the dark matter. Contours labeled with lepton number 
L — 0, L = 0.003, L = 0.01, L = 0.1 are production predic- 
tions for constant comoving density of Qs = 0.24 for L = 0, 
and Sis = 0.3 for non-zero L [2^. Constraints from the CXB 
with the minimal MW halo model are in the solid (blue) re- 
gion, while the maximal MW halo model excludes the ad- 
jacent diagonally hatched region. Also shown are exclusion 
regions from the diffuse X-ray background (green) [g^], from 
XMM-Newton observations of the Coma and Virgo clusters 
(light blue) [68j , observations of Andromeda (M31) in wide 
hatching 40], and limits from the MW by Boyarsky et al. [4l| 
(BMW). The region at rus < 0.4 keV is ruled out by a con- 
servative application of the Tremaine-Gunn bound [l^l. The 
grey region to the right of the L = case is where sterile 
neutrino dark matter is over-produced. The constraint from 
the MW calorimeter soft X-ray background observation is the 
star and arrow, marked "Calor." Also shown is the horizontal 
band of the mass scale consistent with producing a 100-300 
pc core in the Fornax dwarf galaxy [5(. The non-resonant and 
resonant production curves come from Refs. [2^ and [29l ]. re- 
spectively. 



per 2.5 eV bin, so ^ 5 counts within FWHM of 9 eV. For 
Poissonian statistics, for an emission line on top of this 
background, a Str upper hmit is 12 counts (9 counts 
within the FWHM), or 2cr limit of 6 counts. This would 
place a 2a limit between 3 and 11 times the mixing an- 
gle sin^ 29 predicted by the DW model, for the high and 
low mass MW models, respectively, at a sterile neutrino 
particle mass of 2 keV. For the low mass case, this corre- 
sponds to excluding mixing angles at sin^ 26' > 3 x 10^^ 
(see Fig. [3]). At lower energies the detector efficiency and 
the neutrino decay rate both decline, so the constraints 
weaken rapidly. Nevertheless, the above limit for TOs = 2 
keV is similar to our ACIS constraints from a much longer 
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FIG. 4: The upper panel is the spectrum soft X-ray back- 
ground as measured by McCammon et al. [H^]. The lower 
panel shows the atomic line model, detector response function 
(blue) and expected contribution ( x 60) due to sterile neutrino 
decay of the DW production model with the minimal (lower, 
green) and maximal (upper, red) MW halo models. 



observation, which shows the huge benefit of a wide FOV 
combined with a calorimetric energy resolution. 



DISCUSSION AND CONCLUSIONS 



The measurement by Hickox & Markevitch [52, l53| of 
the unresolved cosmic X-ray background, produced by a 
total of 3 Ms observations by the ACIS CCD aboard the 
Chandra X-ray telescope, presents a constraint on the 
presence of a line flux that would be expected due to the 
decay of the sterile neutrino dark matter candidate com- 
prising the dark matter halo of the MW. We find that 
the 2(7 upper limits on the sterile neutrino particle mass 
in the simplest (DW) production mechanism are in the 
range of nig < 2.9 and < 5.7 keV for high and low 
mass estimates of the MW dark matter halo, respectively. 
To be conservative, the higher value for the upper limit 
of 5.7 keV is the robust 95% CL upper Hmit to conclude 
from the CXB analysis in the work presented here. More 
accurate dynamical measures of the MW halo are directly 
relevant for the particle mass constraints presented here. 
In addition, we find that there is a significant limit to 
potential constraints from the Chandra observatory at 
photon energies Ej < 1 keV due to the presence of em- 
mision lines of the low temperature gas in the MW as 
well as the rapidly worsening relative energy resolution 
of the ACIS (as well as XMM EPIC) detector at these 
energies. This limits potential constraints from Chandra 
(or XMM-Newton) on the sterile neutrino particle mass 
to not better than m-s ~ 2 keV. 

The limits presented here from the MW contribution 
to the measured CXB in Chandra are comparable to 
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that from the XMM-Newton observation of Andromeda, 
rus < 3.5 keV (95% CL) for the DW model, by Wat- 
son et al. [Jo] who use a stringent requirement of the 
decay signal to be four times the astrophysical back- 
ground within a bin to place this upper limit. Their 
constraints in the full parameter space (Fig. [3|) were de- 
rived using an energy-averaged flux-to-counts conversion 
instead of the direct spectral fitting, which resulted in 
an approximate power-law constraint. In reality, their 
constraints must weaken at high energies as they do in 
this work, because the XMM-Newton effective area also 
declines sharply with energy. 

The MW constraints by Boyarsky et al. [s^ from the 
XMM-Newton X-ray background data in the full param- 
eter space are comparable to ours. Their more recent 
analysis of the XMM background data i41| resulted in 
3(7 constraints plotted in Fig. 3 (where the horizontal 
gaps correspond to the intervals E = 2.3 — 2.6 keV and 
5.9 — 6.3 keV excluded from their spectral fits). Their 
adopted MW halo model lies between our minimum and 
maximum mass models. While Chandra has an advan- 
tage over XMM-Newton of a more stable instrumental 
background and almost complete removal of the CXB 
in our deep exposures, the solid angle subtended by the 
XMM-Newton FOV is about 10 times greater than ours. 
This results in a proportionally higher expected MW ster- 
ile neutrino signal and more stringent constraints. On 
the other hand, because of the above two reasons, we 
could afford to derive constraints in a slightly wider en- 
ergy range. 

Reimer-Sorenson et 



al. 42| also analyzed Chandra 



CXB observations using the blank-sky data, and derived 
constraints on the MW sterile neutrino flux. Because 
they did not remove the ACIS internal background (that 
dominates at all energies and is well calibrated), their 
constraints on the flux are about two orders of magni- 
tude weaker than ours. 

We have also analyzed a related observation by Mc- 
Cammon et al. [IJ] in the soft X-ray, motivated by the 



high spectral resolution of the observation at low photon 
energies, as well as the large field of view (0.81 ster) pro- 
viding a corresponding two orders of magnitude larger 
MW halo mass in the field of view compared to the 
CDF observations. We find that though this observa- 
tion in itself does not provide competitive constraints in 
the DW production model, longer exposure observations 
with low backgrounds may detect or place constraints 
within a photon energy range {Ey < 1 keV) and par- 
ticle mass range (tos < 2 keV), which is inaccessible 
to current X-ray telescopes. This observation constrains 
rus — 2 keV neutrinos to not have sin^ 26 > 3 x 10"''. 
Although this region is excluded in the DW model due 
to over-production, the constraint may be of interest for 
low reheating temperature or dilution production mod- 
els ^32,, ^3^] . Future observations of the type by McCam- 
mon et al. in the soft X-ray, as well as those possible 
with the Constellation-X observatory at higher photon 
energies [35|, may probe the full parameter space for 
oscillation-based sterile neutrino dark matter production 
models. 
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